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Abstract. We provide an explicit combinatorial expansion for the ground 
state energy of the massless spin-Boson model as a power series in the coupling 
parameter. Our method uses the technique of cluster expansion in construc- 

\^J tive quantum field theory and takes as a starting point the functional integral 

representation and its reduction to an Ising model on the real line with long 

\^0 | range interactions. We prove the analyticity of our expansion and provide an 

£f*l ■ explicit lower bound on the radius of convergence. We do not need multiscale 

nor renormalization group analysis. A connection to the loop-erased random 

l/~) ■ walk is indicated. 
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1 Introduction 

The spin-Boson model is one of the simplest quantum mechanical models 
for the interaction of matter with radiation. Yet, many of its properties, 
especially scattering in the massless case, remain a mathematical mystery. 
It is (formally) given by the Hamiltonian 

H x = (cr z + I)® I + I <8> / w(k) a*(k)a(k) d z k 



Jm.3 47r v / w(/c) 

The Hilbert space is C 2 ® J- where J- is the Bosonic Fock space built 
on f) = L 2 (R 3 ,C,d 3 /c). The creation and annihilation operators, a*(k) and 
a(k) respectively, satisfy the usual canonical commutation relations. The 
matrices a x and a z are standard Pauli matrices. We refer to [35] for a more 
detailed spectral theoretic presentation of the model. For ease of reference, 
we kept the exact same formulation and notations as in the cited article. For 
more background from a physical perspective, the reader may consult [15] 
for a good review of the spin-Boson model as well as its Fermionic cousin: 
the Kondo model. For a review from a mathematical perspective, we rec- 
ommend [39], where many open problems are listed. 

It is well-known that if -4= and £ belong to f), and for real A, then H\ 
is self-adjoint and bounded below. Of particular interest is the massless 
case where u(k) = \k\ and / does not vanish at the origin. Because of the 
infrared singularity at k = 0, the mathematical analysis in that situation 
is more difficult. Once self-adjointness is established, the first step towards 
the understanding of the dynamics of the model is to construct the ground 
state f2cs(A) and its energy E(X) = inf a(H\). As this energy is at the 
bottom of the continuous spectrum when A = 0, the standard perturbation 
theory tools (as one can find say in [II]) do not apply. More sophisticated 
methods are needed for the control of the ground state energy for nonzero 
coupling A. Such a method was introduced by Bach, Frohlich and Sigal in 
their seminal work |1T|. [T8] and it has been further developed over the last 
decade (see, e.g., [14L I15|. fT6] ) . One of its main motivations is to provide rig- 
orous constructive algorithms for the calculation of quantities such as E(X), 
in the case of massless Bosons. Such an algorithm based on the renormal- 
ization group (RG) in operator space, via the Feshbach transformation, is 



in principle powerful enough to show analyticity results for E(X). The full 
delivery on that promise, for the spin-Boson model, has been made only 
very recently in [35] where the construction of E(\) and Qgs W an d the 
proof of their analyticity was obtained. In fact, two constructive algorithms 
are given in [35]. The first one is to follow the RG iterations. The second, 
which is computationally much simpler, can be summarized by saying: put 
an infrared cut-off, compute the Rayleigh-Schrodinger (RS) pertubation se- 
ries then remove the cut-off. At the end of this process, all the terms of 
the series are free of infrared divergences. Moreover, the series converges 
for small A. Unfortunately, the first algorithm is needed to prove these 
statements, which are very far from obvious when looking at the RS series. 
Indeed, by looking at the first few low orders in perturbation, Hasler and 
Herbst found some rather mysterious cancellations of infrared divergences. 
By their RG proof, one knows in an indirect way that these cancellations 
must persist through all orders in A. Prior to |35j, analyticity of the ground 
state energy with respect to A (and other parameters) was established for 
less singular nonrelativistic QED models, using the RG techniques of Bach, 
Frohlich and Sigal, in the article |33j. On Page 580 of the latter, one can 
read: "It seems unlikely that another approach, not based on a renormal- 
ization analysis would yield a result similar to ours". As demonstrated in 
the present article, there is another way ! 

Another expression for E(\) was found in [36J but it is far from explicit, 
nor is it clear (at least to us) how one can use it to prove analyticity. In this 
article, we obtain a completely explicit combinatorial expansion for E(X). It 
is explicit to the point one can see that it is analytic literally just by looking 
at it! By "it" we mean the graphs produced by our expansion as in the figure 
at the end of §4 below. Moreover, the new method we introduce in this article 
does not involve multiscale analysis nor renormalization, despite the massless 
situation. Our method is a direct descendent of the infamous combination of 
the Euclidean strategy plus the cluster expansion technology which worked 
so well in the context of constructive quantum field theory (CQFT). We are 
not the first to use this general approach in the context of nonrelativistic 
QED, see, e.g., [45]. However, we have the advantage of using the latest 
'stable version' of the CQFT cluster expansion software [Tj which comes 
with such optional features as multiscale analysis and renormalization, large 
versus small field analysis, p-particle irreducibility analysis, and could also 
handle microlocal sectorial decomposition as in |48[ 124] . Since we do not 
need these in this article, we keep them in reserve for more dire infrared (or 
ultraviolet) future circumstances. There is also a very promising new set 
of ideas on how to improve CQFT expansion techniques (see [36, 34]), but 



this is still at the 'beta version' stage. These ideas which originated in |50j 
unexpectedly came from the study of renormalization in noncommutative 
quantum field theory (see, e.g., [52]). A good pedagogical entry point in 
the subject of CQFT expansions is [51] (see also the follow up research 
article [55]). 

Our point of departure is Bloch's formula 

E(X)=lim o -llog(n h e- TH ^) 

as well as the Feynman-Kac-Nelson formula followed by integration over the 
Bosons, with the result 

(n,, e - r *n,) = z((A) 2 , T ) (i) 

where 

exp ( - / / X(t)X(s)h(t - s) dt ds 



Z(a,T)=E 



The vector £2|_ is the free vacuum. The function h(s), s 6 1, is defined by 

h{s) = f IZM! e -\sWk) d 3 fc > o . 

J R s uj(k) 

Under the hypotheses on /, it is continuous and satisfies ||/i||l°° = || /=||l 2 - 
The finite time (or volume) partition function Z(a, T) is an expectation over 
a continuous time Markov jump process X(t), t > 0, taking values ±1. The 
boundary condition at t = is free: P(X(0) = 1) = P(X(0) = -1) = ~. The 
wait times between two consecutive spin flips are exponentially distributed 
with unit parameter (i.e., with measure e~*(i£). This is the famous Ising 
model over R with long range ~ h interactions introduced in [11] for the 
Kondo model and in |25] for spin-Boson model. See [59} [26l [58] for a math- 
ematical study of the spin-Boson model via this route. Our main result is 
the following. 

Claim: One has an explicit power series expansion for 

lim — -logZ(a,T) . 
T-»oo 1 

This quantity is analytic, with a lower bound on the radius of convergence 
in a given by 

Rmm = [32 x y/e x max (| |/i| |^oc 



We defer the precise statement of the theorem to §4 since a fair amount 
of combinatorics must be introduced beforehand. Note that we do not yet 
have a lower bound on the radius in A for the ground state energy. This is 
because Bloch's formula requires a nonzero overlap (fi^, £Igs (^)) 7^ which 
is known to hold for small coupling, for instance by the results of |35| . but 
with no explicit range of validity. This will be fixed in [5] where the present 
expansion technique is used for the construction of the ground state itself. 

We will not dwell on the the proof of the reduction ([T|) which is stan- 
dard. One needs to use the Schrodinger representation of T in Q-space (see, 
e.g., [32j[56]) as well as the Feynman-Kac-Nelson formula (see, e.g., [56 | l42 | 
[37] and also |38j for the simultaneous treatment of spin). Note that the 
spin part of the interaction can be made to act as a multiplication operator. 

This is done by a simple conjugation with « = i J , namely, using 

u^ 1 [a z + I)u = I — o~ x for the free Hamilatonian and u~ 1 a x u = o~ z for the 
interaction. The jump process is the one corresponding to the semi-group 

e -t(I-(T x ) _ Fi na riy the vector u^ 1 I I coming from the spin-down choice 

fij,, is responsible for the boundary condition on X(0). 

By contrast, the combinatorics of our expansion will be explained in 
great detail. When used in problems of mathematical analysis, combina- 
torics is a bit like a Toruk [13J -—it can be a powerful friend, but it must be 
approached with respect. The degree of detail serves the purpose of show- 
ing, by example, what we mean by respect for combinatorics. Over and 
over in the course of our proof, we will use relabeling changes of variables, 
continuous or discrete. This is because behind the scene there is Joyal's 
theory of combinatorial species |40[ [j~9] at work. Although we do not re- 
sort to this theory explicitly, in this article, it is the proper mathematical 
framework for much of the combinatorics involved in Feynman/cluster type 
expansions. The two references mentioned only consider applications to the 
field of combinatorial enumeration. The weights used are placeholders such 
as x n where n counts something. In accordance with Sokal's multivariate 
philosophy [57], allowing more general weights such as contractions of ten- 
sors is where the full potential of the theory is. Some steps in this direction 
were taken in [2]. From the point of view of the latter article, the theory 
of species becomes a lingua franca shared by such remote areas as quantum 
field theory and (neo)classical invariant theory. The only difference is that 
of infinite versus finite dimension respectively. More on this can be found 



in gain Ei- 

This article demonstrates that there is a bright future for 'the Euclidean 
strategy + CQFT cluster expansion' method in the area of nonrelativistic 
QED. The first part, i.e., the Euclidean reformulation of the problems posed 
by the interaction of matter with radiation should be straightforward. To 
paraphrase the prophetic words of Paul Federbush [28], the second part is 
where the action is. One can try to use techniques similar to ours, in order 
to treat progressively more difficult and more realistic models. However, we 
think it is more important, at this point in time, to focus on the simplest 
models but try to prove more, i.e., to try to understand scattering com- 
pletely. As advocated in [39] , the massless spin-Boson model is definitely in 
this class. There is another one which is even simpler: the Nelson model 
with quadratic confining potential. We take this opportunity to announce 
the obtention of the exact analogue of the result in the present article, for 
this particular case of Nelson's model (above three dimensions) [8]. We are 
encouraged by how forgiving the combinatorics of these two models are, and 
therefore led to believe that a deeper understanding of scattering is within 
reach. The ultimate goal, if this endevor succeeds, would be to apply the 
lessons learned in such simpler models to the scattering theory of P{4>)2- 
The latter is the missing piece of the glorious CQFT work in the seventies 
(see [12] for a recent update and list of references) . We would like to dream 
that this achievement is possible. 

2 The jump process 

Let N(t), t > 0, be the usual Poisson process with parameter 1. Let B = ±1 
be an independent Bernoulli random variable with F(B = 1) = P(2? = 
—1) = \. We realize the jump process X(t), t > as 

X(t) = B ■ (-1) N ® . 

Clearly, Z(a, T) defined in ([2]) is an entire function of a G C, for any 
given finite time T. Indeed, by the Fubini-Tonelli Theorem one has the 
everywhere convergent power series expansion 



oo 2p 



Y[h(t 2j -t 2j -i) X E[X(*i)...X(t2p)] . 

3=1 



Our first task is to explicitly evaluate the expectation. 

Lemma 1 Let < t\ < ■ ■ ■ < t q be an increasing sequence of times, then 

e -2[|tl-t 2 | + |t3-t4|H \-\tq-l-t q \] jf g j g eyen ^ 







if q is odd. 



E[X(h)---X(t q )] = 

Proof: Write 

X(h) ■ ■ ■ X{t q ) = B" • (-l)"(*i)+-+"(*,) . 
Then, by Abel's summation by parts 

N(h) + - • -+N(t q ) = qN(ti)+(q-l)(N(t 2 )-N(t 1 ))+- ■ ■+(N(t q )-N(t q - 1 )) . 
By the independence of the increments 



E[X(t 1 )---X(t q )]=E[B*\ 



1 



,qN(ti) 



xE 



-1) 



(8-l)(JV(ta)-JV(ti)) 



1 



,tf(t,)-iV(t,_i) 



This vanishes if g is odd because E(l? 9 ) = E(B) = 0. Now assume q is even, 
then 



E[X(ti)..-X(t,)]=E 



-1 



\N(t 2 )-N(ti) 



(-1) 



N(U)-N(t 3 ) 



■••E 



(-1) 



AT(t,)-JV(t,_i) 



Since by the definition of iV, N(ti) — N(ti-±) is distributed according to the 
Poisson distribution with parameter tj — tj_i, 



E 



(_l)M*i)-^(*i-i)l = V( l)fc c -fe-^:-i) & ~ fr- 1 ) 
fc=0 



-2(tj— tj_i) 



and the lemma follows. 



3 Preparation for the Mayer/cluster expansion 

In order to use the previous lemma, we need to break the cube of integration 
over times into (2p)\ simplices. We therefore write 



p=0 y 



2p 
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][[ h (t2j ~ t 2j -i) x Yl e - 2 ^W)-^W-n\ . (3) 

i=i i=i 

We find it more convenient to work with a fixed simplex, so we perform a 
change of variables 

(*<r(l) > ■ ■ ■ > *(r(2p) ) -^ (*1 > • • • > *2p) 

for given fixed permutation <r. We get 

2p 



Z{a,T) = Y(°Y± Y [ T\dt 



[I h(t a - 1{2j) - i (T - 1(2j _ 1) ) x n e-^-^l . 
i=i i=i 

We will also change variables c -1 — > a in the sum over permutations in 
order to lighten the notations, i.e., have t a ( 2 j) ~ t ff (2j-i) instead of t ff -if 2 j) — 
*o— 1 (2j-i) as arguments of the /i functions. 

Now let us recall some elementary notions of graph theory. We will use 
the notation [n] for the set of integers {1,2, ... , n}. When considering a 
simple graph G = (V, E) with fixed vertex set V, we will identify G with 
its edge set E seen as a set of unordered pairs in V. Namely, let V(V) 
denote the power set of V, and let V (2) = {A £ V(V)\ \A\ = 2} where | • | 
denotes cardinality. We will view E (as well as G) as a subset of V^ 2 ' . If 
G is a multigraph with (unlabeled) repeated edges but no loops ('tadpoles' 
for physicists), we will instead view G ~ E as a multiset of edges, i.e., an 
element of N v . The degree of a vertex is the number of edges incident 
to v and it is denoted by deg G (w). A A;- factor on V is a graph G such that 
deg G (v) = k for all v G V. A perfect matching is a 1-factor. 

Now let us rewrite Z(a,T) in terms of perfect matchings. For p > and 
o" G ©2p, let -P(a) be the perfect matching on [2p] given by 

P(a) = {{a(l), a(2)}, {a(3), a(4)}, . . . , {a(2p - 1), a(2p)}} . 

The one corresponding to the identity permutation a = Id will be denoted 
by P* which is an analogue of a base point. Now by collecting permutations 
a which give the same perfect matching P we have 



oo „ 2p 

n D J0<ti<-<t2„<T ;_ 1 



2p 

p= Q ;■ J0<t 1 <-<t 2 p<T\ 



n h(t a - tb )x n *- 2ita - tbi 

{a,b}eP {a,b}eP* 

where the second sum is over all the ^j perfect matchings on [2p\. Indeed, 
for any P there are 2 p p\ permutations a for which P(cr) = P. 

Now we will partially undo what we just did. By writing the number 1 
in a more complicated way, we have 

Z(a,T) = f>Wi E l J [ ft* 

p=0 P \ P - r G 6(P.) J J 0<ti<-<t2 P <T i=1 

n h(t a - tb )x n e- 2ita ~ tbi 

{a,b}eP {a,b}&P, 

where ©(P*) denotes the set of permutations of the set P*. Now by Fubini's 
Theorem 



2p 

(it. 



2(*r)=E^r e / n 

p= o p " Te6 (p t ) ,/o<tl< '" <t2p<T *=i 

Ie n h(ta-tAx n e 

{ P {a,b}eP J {a,b}GP t 



■2\ta-t b \ 



Given r £ ©(P*) define an induced permutation <t[t] £ ([2p]) by 

a[r](2i-l) = minr ({2i - l,2i}) 
cr[r](2i) = maxr({2i-l,2i}) 

i.e., this is the permutation on the underlying set of points [2p] obtained by 
rigidly permuting the pair-blocks of P* according to r. Clearly, P*, seen as 
a set partition of [2p], is invariant by a[r] so 

TT g— 2|t„— 1 6 | _ TT g- 2 l*<T[T]-l(o)- t <7[r]-l(6)l 

{a,6}GP» {a,6}eP» 



Besides, 



e n ^ ~ ^) = e n ^(^M-Ha) - ^M- 1 ^)) 

P {a,b}eP P {a,b}eP 



since one can sum over the perfect matching P T = {a\r\(A)\A G P} instead 
of P. This is just another (discrete) change of variables P — > P T due to 
relabeling. Now do yet another change of variables 

which results in the equality 



/ n* Ie n h(ta-tAx n 

J0<ti<-<t 2p <T i=1 ^ p {ab}£P J { a ,b}eP, 



e -2\t a -t b \ 



2p 

-2\t a -t b \ 



= [ n* e n m*.-*6) >< n e ~ 2|ta 

J0<t CT[r](1) <-<t CT[r](2p) <T - =1 ^ p { aj 6} eP J { a , b } e p„ 

This way we keep the integrand the same for all r's. The sum over these 
permutations simply rebuilds a larger domain of integration. Namely, 



oo p , 2p p 



2p 



{the intervals "If \ 

[%-i>%] >me n h(t a ~t b )\x n e -2i*«-* & i 
are disjoint J [ p {a,b}eP J {a,6}eP, 

where !{■■■} means the sharp characteristic function of the condition be- 
tween braces. As a result 

oo p ,. 2p p 

z ^ t ) = e fr e / n T 12 n «* n ^--i < m 

p =o p - p J l°> T r p i=i j=i 

x n Kta-t b )x n c- a i*--*»i 

{a,6}eP {fl,6}eP, 

X J] l{[*2t-l,*2i] n[t 2 j-l,*2j]=0} (4) 

l<i<j<p 

i.e., we are ready for a Mayer/cluster expansion. 



10 



4 The expansion 

Firstly, given a perfect matching P and the fixed matching P*, build the 
2-factor P* + P obtained by superimposing the edges of P with those of P* . 
Note that P* + P can be a multigraph since some edges can become doubled. 
Clearly, P* + P breaks up as a disjoint collection of cycles of even length. 
An example is 

10 - ID 




13 ^— . 







where we represented the edges of P* by thick lines and those of P by thinner 
ones. Let 7r(P* +P) be the partition of the vertex set [2p] made of connected 
components of P* + P, i.e., the collection of supports of these cycles. Since 
P and P, also are set partitions, another way of seeing 7r(P if + P) is a P* VP, 
the join in the partition lattice of [2p] where the ordering ^ means "is finer 
than or equal to" . We also denote by n(P* + P) the partition of P* induced 
by the partition 7r(P* + P) of [2p]. Namely, 

n(p* + p) = {x e v{P*)\ u A6 x a g tt(p, + p)} . 

This is the same as the connected component partition of the graph on 
P* obtained from P* + P by contracting the pairs of P* into points (and 
removing loops if any). For a pair A = {2i — l,2i} £ P* we will use the 
notation tA for the interval [t2i-ij^2i]- 
Now write 

I] I{[*ai-i,t2i]n[t^_ 1 ,t ! y] = 0}= H i{t A ntB = 0} 

l<Ki<p {A,B}ePi 2) 
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where the function 



(2) 



F . R n(p t +p) 

s = ( S {x,y}){x,y}en(p»+P)(2) l— ► F ( s ) 

is defined as follows, and 1 is the parameter vector with all components 
equal to 1. By definition, we set 



f(p)= n n i{tAnt B =0} 

xen(p,+P) \{A,B}exW 



A&X,BeY 



* n 

{x,y}en(p»+p)( 2 ) 

Now use the Brydges-Kennedy-Abdesselam-Rivasseau formula [2IJ |9] 
(see [1] for a gentle introduction and detailed proof of this identity). The 
outcome is 

on n(P„+P) 

Here the sum is over all forests J on the vertex set n(P* + P), i.e., simple 
graphs with no cycles. For each edge I £ $ one associates a real variable 
ui integrated between and 1. The integrand is the partial derivative of 
F with respect to the entries S{x,Y} such that {X, Y} £ $, evaluated at 
the vector s(£,u) = (*>(& «){x,y}) {x,y} 6 n(P*+P)( 2 ) defmed b ^ the following 
rule. 

• If X and Y belong to different connected components of 5, then 

• Otherwise, if X and Y are in the same (tree) component of #, then 
s($, u){x,Y} = nun; Ul wn ere I ranges over the unique simple path 
joining X and Y, in the forest $. 

The outcome of the formula involves, for each {X, Y} G J, derivatives 
of the form 



9s { X > Y } ArzX ^_y 



II (l-s { x,y}l{iAnt B ^0}) 
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= - Yl HtAnt B ^(/)}x Yl (l-»{x,y}l{U'nt B '^0}) ■ 

AeX,BeY A'ex,B'ev 

(A',B')^(A,B) 

Namely, for each 'macro' edge {X, Y}, one needs to choose an underly- 
ing 'micro' edge realization {A, B}, and sum over these choices. A better 

(2) 

rewriting of the outcome is in terms of the collection F C P, of such edges 
{A,B}. The result is 

n i{[*m-i, t 2i \ n [t 2 j-i,t 2j ] = 0} = yi I n dvi 

l<i<j<p F •'[°' 1 ] F ZgF 

H (-l{t A C\t B ^%})x \{ {l-r(F,v) {A>B} t{t A nt B ^®}) • 

i^B}€F {A,B}ePi 2) \F 

(5) 
The sum is over simple graphs F on the vertex set P* which can have at most 
one edge {A, B} with A £ X, B £ Y, for any given pair {X, Y} of blocks 
of the partition T1(P* + P). No edge {A, B} is allowed if it is internal to 
some block of the partition II(P* + P). Finally, we require that the induced 
simple graph F on T1(P* + P) be a forest. The graph F is made of all edges 
{X,Y} £ II(P* + P)( 2 ) for which there exists {A,B} £ F with A £ X and 
B £ Y. If I = {A, B} £ F we use the notation I for the edge it thus induces 
in F. For each I £ F one associates an integration variable v\ £ [0, 1]. Now 
r(F, v) is the vector 

( r ( f ^){A,B})) {AB}eP m 

defined as follows. 

• If A, B belong to the same block X of n(P*+P), then r(F,v)r AB \ = 1- 

• If A, B belong to different blocks X, Y of II (P* + P) respectively, 
then r(F, v) = s(F, u){x,Y} where u is the parameter vector defined by 
u? = v\ for all I £ F. 

We can now incorporate identity (|5|) into our last formula (HJ) for the 
partition function: 

oo p „ .. 2p p 

z ^- t ) = e^ee / l]P n*> L _,n* n n*-> < *«} 

x II (-l{t A nt B /0})x [] (l~r(F,t;) {AB} l{i A nt B ^0}) 

{AB}GF {A,B}eP< 2) \F 
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x n e - 2 i"ix n h( ta -t b ) 

{a,b}eP t {a,b}eP 



(6) 



The integrand/summand completely factorizes according to the connected 
components of the multigraph P« + F on P*. The multigraph P* is made 

(2) 

of the edges {^4, B} 6 P* such that there exists a £ A and b <E B with 
{a, 6} £ P. If two such pairs exist (no more can occur) then the edge {A, B} 
is doubled. Finally, P^ + F is the superposition of the edges from P* and F. 
For instance, after contracting the pairs, a possible configuration is 



9,10 



1.2 



13,14 



15,16 




11,12 



17,18 



where P" is represented by thin lines and F by squiggly lines. 

Let nn(P" + F) be the partition of P* made of connected components of 
P" + F. The crucial property of the identity used is that if A and B are not 
connected by P" + F, then r(F,v)^B} = 0- Then complete factorization 
means, in a first step, 

oo p 

^«> T ) = E^EE II m*>x) 



p=0 



p 



p f xenn(P»+F) 



where 



A(P,F,X) = f T7 dvi [ If du 



leFx 



]J HtminA < Wa} X ]J (-l{t A H t B 7^ 0}) 

Aex {A,B}eF x 

{A,B}eX( 2 )\F x 

x n e-^-^x n ^-t 6 ). 

{0,6)61 {a,b}&P x 
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Here, Fx = F n X^ 2 > is the restriction of F to the component X. Likewise, 
Px = P* H X^ 2 > (as a multiset, i.e., the restriction of the corresponding 
multiindex in N p * to X^ 2 >). It is easy to see that the computation of 
V (^, V ){A,B} is purely local to the component X containing the pair {A,B}. 
Namely, it only involves the knowledge of the data Fx, Px, and vx = 

(vi)ieF x - 

In a second step, one also has the factorization 

ee n ap,^x)= n B(x) 

p f xenn(ptt+F) xenn(ptt+F) 

where 

£(*) = EE^( P *> F *> X ) 

Px fx 

since the combinatorial constraints on Px, fx are local to X, i.e., do not 
depend on what happens outside X. 

Finally, B(X) only depends on the cardinality of X, i.e., B(X) = C\x\ 
where, for any p > 1, C p is defined as follows. We let 



2p 
P F connecting J [°' 1 l F leF J [0,T] 2 P i=l j =] 



_ p p 2p p 

c p=n e / U dv u n ^ n ^j-i < %> 



x H (-t{t A nt B ^Q)})x H (l-r(F,v) {A)B} t{t A nt B ^®}) 

{A,B}eF {A,B}eP< 2) \F 

{a,b}£P* {a,b}eP 

where the constituents of this formula are defined exactly as for (|6|) except 
for one difference: now we have the extra condition on F which is that it is 
connecting, i.e., IHI(P" + F) = {-P*}. All we just did is use a bijection of X 
with the set P* of consecutive pairs of [2p] with a new p given by p = \X\, 
in order to relabel everything. 
Now, 



oo 



^ T ) = E^rE E c \x,\---c\x kl 

p=0 y ' k>0 {x lt ...,x k } 

partition of P* 

Ea p y^ 1 v^ J the Xi are disjoint 1 

-pll^kl 1. 1 { andUl 1 = P» }C\x,V-C\ Xk \ 

p=0 r k>0 (X 1 ,...,X k )eV(P*) k k J 

15 



p=0 ^ fc>0 Pi,.-,Pfc>l U=l ) y y 

by the multinomial theorem. 

Up to now this was pure algebra and there was no issue of convergence. 
Now some analysis is needed in order to write 



*^=Ee e e« p 4e 



fc>0 pi,.--,Pfc>l p>o 1. 1=1 J 

-Ei E (^)-(^)- P (E^). 

fc>0 Pi,-,Pfc>l v y v y \p>i / 

By the the Fubini-Tonelli theorem, all we need in order to justify these last 
steps is the absolute convergence estimate 

^ \a\P\Cp\ 

/ i — < °° 

p>l L 

which will follow from ([7]) below. This granted, we have 

log Z(a,T) = ^^ E - 
/ -^ pi 

Now use the fact that 

2p 



'[o,n-'<=i 

with an integrand 



C p = \\dti C p (ti, . . . , *2p) 

J\0,T] 2 p 7LT 



c P (ti, . . . , t 2p ) =j2 yi / F n dv i n ^i^-i < * 2 j 

P F connecting ^I ' 1 ] 1 " «eF j=l 



} 



x II (-i{t A n t B / 0}) x [] (i-r(F, w ) M)B} i{i A nt B /0}) 

{A,B}6F {A,B}GPi 2) \F 

x n e-^-^x n h(t a -t b ) 

{a,b}eP, {a,b}eP 
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which is defined for any i's in R and is translation invariant: 
C p (h + r, ... ,t 2p + T) = C(t ll ...,t 2p ), Vr. 

[rp rj-i -i 2iD 
— 2", g- J , as 

well as impose a 'pin' t\ = in the argument of C p . Hence, 

2p 



logZ(a,T)= V^r / , TT^C p (0,i 2 -ii,...,t2p-ti 

p>l ^' ■ / | _ T>2j i=l 
y^—r dti fTdsj C p (0, s 2 ,... ,s 2p ) 



p>l ^ "-^2 "^ 1 = 2 

2p 
X 
i=2 

after changing variables 



(t 2 — h, . . . , t 2p — t±) — > (s 2 , . . . , s 2p ) 



and implementing the domain specifications by sharp characteristic func- 

T 



tions. Now change variable to w = ^ to get 



log Z(q,T) _ ^«p H '■ 2/ ' 



V]— - / ck; / TTdsj Cp(0,s 2 ,---,S2p) 



T 

p>l L " 2 "^ 1=2 

xfji|-rQ + «;) < Si <rQ- w 

For almost every u/, the characteristic function goes to 1 pointwise, as T 
oo. By the Lebesgue Dominated Convergence Theorem, we will have 



log Z(a,T) v-^ " p I rr , ^ /„ » 

hm — i ; = -E~r / n^ c p°' f2 '---' t2 p 



as soon as one can prove the key estimate 

2p 



r = V^-/ TT^ |Cp(0,i 2 ,...,i 2p )|<oo (7) 



p>\ L "^ i=2 
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which will be done in §5. 

One easily sees from the previous manipulations that the estimate 
implies the earlier one needed for finite T. Namely one has 

E|a| p |CJ ^ „ 
^i^<T-r . 

v>i Pl 

At last, we can state our main theorem with full precision. 

Theorem 1 We have 



a. - loeZ(a - T) " ' ' * 

p>i 



e^e e t.n^t^n* 

p>l * P F connecting J I ' 1 ! ZeF J K 8=2 



[J lfe-i < fej} x [] (-l{unt B ^0}) 

3=1 {A,B}eF 

x [] (l-r(F, V ) {AjB} i{t A nt B /0}) 

{A,B}ePi 2) \F 

{a,b}£P t {a,b}eP 

where t\ = by definition. The limit exists and the sum over p converges 
as soon as \a\ < R m { n with 

R min =[32x\/ex max(||/i||z,°o,||/i|| L i)] _1 

which is a lower bound on the analyticity radius. 

An example of graph produced by the expansion (in the contracted pairs 
representation) is 

9,10 



1.2 



13,14 

11,12 



15,16 

17,18 

where we did not include the interpolated hardcore conditions for a clearer 
picture. 
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5 The key estimate 

We now prove the previous key estimate usually called a polymer bound in 
the CQFT literature, or the Kotecky-Preiss criterion. We start from the 
trivial bound 

ro<£f£ e / 0I]F n^t_,n* 

p>l y P F connecting J l U,i l ieF JK i=2 

P 

Yl l{t 2j -i < t 2j } x Y\ t{t A n t B / 0} 

i=l {A,B}£F 

x 11 \l-r{F,v) {AB} l{t A r\t B ^®}\ 

{A,B}ePi 2) \F 

{a,fe}eP» {a,6}eP 

Since, by definition, the interpolated couplings r(F,v)i A m belong to the 
interval [0, 1] we can use 

\l-r(F,v) {AB} l{t A nt B ^(b}\<l (8) 

and also drop the integral over the v's. Thus, 



r »sE^rE E iW.F) 

P>1 P F connecting 

where 

v( P , p, f) = / rr dti n i{t 2j _! < t 2i } x rr i{t a n t B + 0} 



x n e - 2|ta - tb| x n ^a-h). 

{a,6}eP, {a,b}eP 

o 

For each cycle of the multigraph P* + P, delete one edge of P. Let P be the 

o 

graph P with these edges removed. Clearly, | P \ = \P\ — \ir(P* + P)\. Note 
that we are simply choosing an opening of all the cycles. No summation 
over these choices is made since we are doing an L°° rather than an L 1 
bound. Now it follows from all the previous combinatorial definitions that 
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T = (py + F is a spanning tree which connects the vertex set P*. Here, (p)" 

H ° H ° 

is defined in the same way we did for P ff earlier. Also note that | (P)' | = | P | • 
We now have 

r- 2p p 

v(p,p,f) < \\h\fiL p - +n x / fr^ n%4<t 2J } 

J^-^t= 2 1=1 

x 11 l{t A n t B + 0} x 11 e -2|*»-*bl x J] fr(t a - t 6 ) 

{A,B} 6 F {a,6} 6 P, {o6}ej ° 

and recall that t\ = by definition. We now simply use an inductive L l -L°° 
bound or 'pin and sum' argument for the integration over the i's, by following 
the tree. These integrations have to be done in pairs t^i-x-, t^i- Indeed, what 
follows should be thought of as a 'sum' over all possible intervals t^, for 
A G P*. The element {1,2} e P* is chosen as the root of T and one 
bounds the integrals over the t^'s starting from the leaves, conditioning on 
their ancestors, and then progressing towards the root. The procedure is 
somewhat subtle since it is a mixture of straightforward pin and sum with 
L 1 decay and a generalization of Cammarota's argument [22] to polymers in 
the continuum instead of on the lattice. A possible choice of cycle opening 
for the example in the last picture is 



9,10 




11,12 



15,16 

k 17,18 

where we indicated the orientation towards the root 1,2 for the edges in the 
remaining tree T. 

Let us first explain the integration over £2^-1 and tii when A = {2i— 1, 2i} 
is a leaf. 
1st case (e.g., pair 17,18 of example): If the edge joining A to its 

ancestor B in the rooted tree T comes from (P)", then one needs to bound 
an expression of the form 

f dt 2i -i dt 2i i{t 2i -i < t 2l } e - 2 i^-i-<2»i m a _ t6 ) 
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where a £ A and b G B. In this case, first integrate on the t labeled by the 
other element of A, then integrate on t a . The result is ^H/iIIl 1 - 
2nd case (e.g., pair 5,6): If the previous edge conies from F instead, then 
one needs to bound 

/ dt 2i -i dt 2i Ht 2i -i < t 2i } e - 2 i^-i-^i i{t A n t B + 0} . 

Jm. 2 
Let B = {2j — 1, 2j} for some j ^ i- Performing the change of variable 

J £ = Tjfe -hi-i) 

1 V = -^(t2i + t2i-l—t2j — t2j-l) 

the previous integral becomes 

^ / eft> e' 2 ^ = / ^ 2(e + 1 {t2j _ i2 ._ l)} e -2^ 



o J-Z-^(t 2j -t 2j -i) Jo \/2 

= 1 (1 + 2(t 2i - t 2j -i)) ■ 

Therefore, when we later come to doing the integration over ts we have an 
extra factor 1 + 2(t2j — ^2j-i) to control per offspring born from an F edge. 

Let us now consider the integral over t A where A is a vertex deeper in 
the tree but distinct from the root. Again one has two cases to consider. 
1st case (e.g., pair 7,8): Same as previous 1st case but with q offsprings 

connected to A by (P)* edges. One has to bound 



/ dt 2i -i dt 2i t{t 2i -i < t 2i } e - 2 l*2-i-^l h(t a - t b ) (1 + 2(t 2i - t 2 , 

JR 2 



-l)) q 



We introduce a parameter 7 > to be fixed later and we write the trivial 
inequality 

(1 + 2(t 2i - t 2i -i)) q < 4 • e^ 1+2 ^- t2 ^» . 
Now, provided 7 < 1, we can proceed as before and get a bound 

^ f dt 2i -l dt 2i t{t2i-l < t 2i } e -2(l-7)|fai-i-*«l h{ t a - tb ) 

l q Jr 2 

~ 2(1-7) 
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2nd case (e.g., pair 13,14): Same as previous second case but with q 
offsprings connected to A by F edges. One needs to control 

/ dhi-i dt 2i l{t2i-i < t 2i } e^**- 1 "* 2 *! l{t A n t B + 0} 

Jm. 2 

x{l + 2{t 2l -t 2i ~i)) q 
= / dU dn e- 2 ^(l + 2V20 q 



V2 
"OO 



poo _ 1 

= J di e- 2 ^ (1 + 2V20 9 x 2(i + -j=(t 2j - t 2j -i)) 
d£ e~ 2 ^ 7 -«g! e^+ 2 ^ x 2(C + -Lfoy - *y-i)) 



g7^y <?n' 



4(l- 7 )2 
Using 7 > 0, we get an upper bound 

q\^f~ q e y 



-[l + 2(l- 1 )(t 2j -t 2j - 1 )} 



4(1 -7) 2 



(1 + 2{t 2j - t 2j -i)) 



This is the crux of our proof: one reproduces a factor 1 + 2(t 2 j — t 2 j-±) of 
the same form as that coming from the offsprings. Therefore, the induction 
works. 

Finally, for the root A = {1,2}, there is only one case to consider. Let 
again q be the number of offsprings attached to it by F edges. Given that 
t\ = 0, the quantity to bound is 

f dt 2 1{0 < t 2 } e~ 2 ^\ (1 + 2(t 2i - t^)) 11 < |)f-^ . 

Jr 2 (l-7) 

As a result of this discussion, we have 

2?(p,P,F)<||/ l ||] L l P * +p)l x||/ l |gf ),| x J] [^T^W] 

AeP* 

where q& denotes the number of F offsprings of vertex A and where 
5(7) = max 



2'4'2(1-7)'4(1- 7 ) 2 '2(1-7)J " 
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The redundancy recapitulates the previous list of cases. In fact, 



6(y) 



when < 7 < 1. Since 



4(1-7)5 



x max{l,2(l — 7)} 



J2 QA = \F\ < |T| <p 



AeP, 



and 



we have 



|tt(P* + P)| + |P| = |P| 



P(p,P,F)<[max(||/i|| LO o,||/i|| L1 )] |P| x<5(7) |P * l X7- p x \{ q A \ 

AeP* 

Since \P\ = \P*\ = p, we therefore have 



r ^Ei 



p>\ 



\a\6(j) 

7 



max 



e e n 



<7A! • 



P F connecting AgP* 

Now recall that the qA depend on p, P and F. However, it also depends 

o 

on the noncanonical choice of cycle opening P, but only via the resulting 
tree T. Therefore a more precise writing of the right-hand side of the last 
equation involves 



mm 



E E m i 

P F connecting 



II <M(P,F,T)! 

L4eP, 



sE E E> 

P F connecting T 



P,F,T 

are compatible 



x 



deg T ({l,2})! [] (deg T (A) - 1)! 



AeP* 
A#{1,2} 



where compatibility means that T can be obtained by some choice of cycle 
opening. Note the use of F C T to bound the g's in terms of the vertex 
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degrees in T. Note that the root {1,2} G P* is special in this regard. By 
Fubini, the bound becomes 

£deg T ({l,2»! JJ <d^(A)-l).x£l{ £^. 

T AGP, pf *- 



A#{1,2} 



<4^deg T ({l,2})! J] (deg T (A)-l)!. 

T Aep, 

A^{1,2} 

Indeed, counting the pairs (P, F) which are compatible with a fixed tree T 
can be done by summing over the number k = |7r(P* + P)|- This determines 
the cardinality of F C T, i.e., |F| = k — 1. Finding F accounts for a factor 

of I i ) • Then (p)' = T\F is determined. Since the connected com- 

O n 

ponents of (P)" are open cycles, there is only one way to close them back, 
so P* is known. However, finding P costs an extra 2 P due to the choice of 
micro realizations {a, 6} £ [2p]( 2 ' of the macro edges {^4, B} of P*. One gets 
a 2 P bound instead of 4? because P is a perfect matching. In sum, 

^—* \ are compatible J ~~ *-^ \ k — 1 J 



At last, we can write 



y- 






max(||/i|| L oo,||/j|| L i) 



1 



X 



4|a|(5(7) 

max(||/t|| Loo ,||ft|| L ^ 

^ P>2 ^ 



£deg T ({l,2})! J] (deg T (A) - 1)! 



T A6P, 

A#{1,2} 



4|«|<5(7) , 

max(||ft,|| LO o,||/i|| L i, 



x 



(p-2)! 



E 7^-T&7^W x ^ 2 -l)!---K-l)! 



^.tr^^- 1 )'---^- 1 )' 

£d i= 2p-2 



by Cayley's theorem [601 Theorem 5.3.4] which counts connecting trees on 
a labeled set with fixed vertex degrees di = deg T ({2z — 1, 2i}). 
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The last sum is bounded using 

£ (p-2)!d 1 <(p-l)! 2 1 



d-^ ,. . . ,dp>l d^ , . .. ,dp>l 

12d i =2p — 2 I2d i =2p-2 



(/,-!)!( ^ ) £(„-l)M". 



Finally 



with 



r <E w 



P >i p 



K = 32xmax.(\\h\\ L <x>,\\h\\ L i)x inf 

0<7<1 



5(7) 

7 



<32x Vexmax(p||Loo,||fc|| L i) 

if one takes 7=2 f° r simplicity. 

As an illustration of the quality of these bounds, let us pick f(k) 
H{|^| < 1} and a; (A;) = |&|. In this case 

,, / 1 Z" 1 

50 = — — \\ L 2 = 4-7T / kdk = 2tt 
V^ Jo 



and 

/'OO /»1 

||/i|| L i =8vr / ds dkk e - sk = 8ir 
Jo Jo 

and therefore we have a lower bound on the radius of analyticity in a 

(256vr v / ^)" 1 ~ 7.54 x 10" 4 . 

Since a = (j^) , this translates (provided Bloch's formula holds) into a 
lower bound for the A radius of analyticity ~ 0.34. 

6 Closing remarks 

Remark 1: One could presumably improve the estimates, in the spirit 
of [30], using the remaining hardcore constraints which were thrown away 
in (J8j) . In this case one should use the Rooted Taylor Forest Formula of [9] 
which forces offspring exclusion. Then one should follow the presentation 
in [27]. The rooted forest formula has more of a Fermionic flavor. It is 
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related to the matrix-tree theorem and is closer to Penrose's lemma for the 
convergence of the Mayer series |J7] . 

Remark 2: Surprisingly, in two dimensions, the obtained series still makes 
sense, order-by-order, even without infrared cut-off. Its Borel transform 
has a nonzero radius of analyticity. This can be shown by a more refined 
analysis along the lines of the previous proof. In this case h(s) ~ - for large 
s. Therefore, one cannot open the cycles. The integral over the t's belonging 
to such a cycle of length k (the number of P* points involved) converges but 
produces a k\ (see, e.g., [39J §3]). In two dimensions, mfa(H\) = — oo and 
we do not know what the term-by-term infrared-finite perturbation series 
means in that case. 

Remark 3: There are similarities between our expansion and the one found, 
e.g., in [54] . and for a good reason: both essentially are resummations of the 
same underlying Feynman diagram expansion. Indeed, the massless spin- 
Boson model is a typical example for which the results of |54j can be applied. 
There are also important differences. The article |54] uses a real-time set- 
ting whereas our article is framed in the Euclidean setting. Also, the cited 
article uses a decomposition into cubes or rather intervals as is traditional 
in Bosonic CQFT. Our approach does not use such a decomposition and is 
closer to the Feynman diagram perturbation series. We decouple abstract 
vertices of such diagrams rather than cells in a decomposition of space- 
time. A similar distinction can be made between existing techniques for the 
treatment of Fermionic theories. For instance, early cluster expansions for 
Fermions used decompositions into cells (see, e.g., [31], EH])- Later, expan- 
sions which decouple Feynman diagram vertices [441 [10] allowed to achieve 
the same convergence bounds in a simpler way. 

Remark 4: The combinatorial structures which feature in our expansion 
for the spin-Boson model are more than strikingly similar to the "canvases" 
of [6] which are related to the loop-erased random walk (LERW). We suspect 
a model mapping such as the one in [20J is lurking behind this similarity. 
Here the LERW seems to be one on the real line where a particle alternately 
'marks its territory' with an interval of exponentially distributed length, 
and diffuses according to the function h. The self-avoidance (or rather loop- 
erasure prescription) is in terms of this previously marked territory. We 
believe a similar LERW model on R 2 instead of R could help bridge the gap 
between the Schramm-Lowner Evolution theory [55] and conformal field 
theory in its Coulomb gas formulation (see, e.g., |23j). 
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